Master equation approach to optomechanics with arbitrary dielectrics 
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We present a master equation describing the interaction of light with dielectric objects of arbitrary 
sizes and shapes. The quantum motion of the object, the quantum nature of light, as well as 
scattering processes to all orders in perturbation theory are taken into account. This formalism 
extends the standard master equation approach to the case where interactions among different 
modes of the environment are considered. It yields a genuine quantum description, including a 
renormalization of the couplings and decoherence terms. We apply this approach to analyze cavity 
cooling of the center-of-mass mode of large spheres. Furthermore, we derive an expression for the 
steady-state phonon numbers without relying on resolved-sideband or bad-cavity approximations. 



I. INTRODUCTION 



In quantum optomechanics, light is used to cool and 
control the mechanical motion of massive objects in the 
quantum regime [TU!] . These systems have potential ap- 
plications in quantum information 5-8 , metrology [9]- 
112] , and can even be used in experiments testing the foun- 
dations of quantum mechanics [T5HT5] . In the broad re- 
search area of cavity quantum optomechanics two classes 
of systems can be distinguished: the reflective case, re- 
alized in deformable Fabry-Perot resonators [THl [T7] or 
microtoroidal cavities |18j . and the dispersive case, like 
in the membrane-in-the-middle configuration [19H21] or 
in optically levitating nano-dielectrics [22TI2B]. In the lat- 
ter, the dimension of the object along the cavity axis 
(i.e., the width of the membrane or the diameter of 
the nanosphere) is typically much smaller than the op- 
tical wavelength. This implies that the dielectric can 
be treated as a dipole with some induced polarizabil- 
ity [23J 121] • The problem is akin to that of single 
point particles, like atoms or ions, in the weak excita- 
tion regime. Thus, the theory and methods that have 
been developed in the context of laser cooling, trapping, 
and manipulation of single atoms and ions can be directly 
applied to optomechanical systems (see, e.g., [27H29] for 
some expository articles). In particular, sideband-cooling 
techniques [30J [3T] have been successfully employed to 
achieve the ground state in a nano-optomechanical sys- 
tem [3J[33] (see also [3J[35]). 

The control that is being achieved in dispersive quan- 
tum optomechanics opens up the challenge to explore the 
physics of larger objects. While this is certainly within 
experimental reach |36j . the existing quantum theories 
are not applicable since the dielectric object can no longer 
be considered as a simple dipole. In contrast, for sizes 
comparable or larger than the optical wavelength, multi- 
scattering processes within the dielectric have to be taken 
into account. As it is well-known from classical nano- 
photonics, they give rise to a modification of the forces 
experienced by the system, as well as other interesting 
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phenomena [37]. 

In this article, we present a quantum theory describing 
the interaction of light with the center-of-mass of non- 
absorbing dielectrics of arbitrary shapes and sizes. In 
particular, we derive a master equation for the motion of 
the particle and the cavity mode. This method consid- 
ers the full scattering process by linking the coefficients 
of the master equation to the scattering matrix. This al- 
lows one to use the tools and techniques developed in the 
context of classical nano-photonics to determine the evo- 
lution of the quantum system. These include advanced 
numerical techniques, like the discrete-dipole approxima- 
tion [35], the T-matrix method [3S], or, for some special 
geometrical shapes, even analytical solutions, like the so- 
called Mie solution 



We apply this approach to a dielectric object in a 
high-finesse optical cavity and obtain an equation that 
is quadratic in the field operators of the center-of-mass 
and the cavity mode. For this general master equation, 
the steady-state phonon numbers are derived without 
relying on the common resolved-sideband or bad-cavity 
approximations. We illustrate this solution by analyz- 
ing the problem of laser-cooling in cavity optomechanics 
with optically levitating dielectric spheres of diameters 
comparable or even larger than the cavity mode wave- 
length. While ground-state cooling can be achieved for 
spheres much smaller than the wavelength, the minimal 
phonon numbers attainable for larger spheres oscillate 
around values of n m j n rs 500. Note that this approach as- 
sumes that all photons are scattered into the bath modes, 
which is justified for objects that are not adapted to the 
cavity geometry, such as spheres or cylinders, but not for 
membranes 1431. 



The manuscript is organized as follows: In Sec. [TTJ we 
describe the system, listing the assumptions and defining 
the Hamiltonian. Following this, we present the main re- 
sult of this manuscript in Sec. |III| a master equation de- 
scribing the interaction between light and the motion of 
arbitrary dielectric objects. First, the effect of the pres- 
ence of a dielectric on a free electromagnetic field is dis- 



cussed in Sec. Ill A where the corresponding scattering 
equations are solved. Based on this, we derive a general 
master equation describing the joint dynamics of the cav- 
ity mode and the center-of-mass motion of a dielectric in 
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Sec. Ill B| The description of the optomechanical setup 



is obtained by assuming the Lamb-Dicke regime and a 
strong driving field in Sec. |III Cj In Sec. |IV| we investigate 
the possibility to cool dielectric objects deriving a general 
theory beyond the common approximations to determine 
phonon numbers in the steady state. Finally the optome- 
chanical parameters for levitating spheres are determined 
in Sec.|Vj In Sec. [VI] we draw the conclusions and give an 
outlook on further directions. In App. [A] we review the 
classical solution of the scattering equations. Following 
this, in App. [B] we describe how the full spectrum of the 
electromagnetic field after interaction with a dielectric 
can be determined within an extended Wigner-Weisskopf 
approach. In App. [C] we supplement the derivation of the 
master equation in the optomechanical case described in 



Throughout the first part of the paper 
Sec 



III 



Sec. Ill C The manuscript is rounded off in App. O by a 



discussion of the small-particle limit. 



II. PHYSICAL MODEL AND HAMILTONIAN 

In this section we describe the system consisting of a 
dielectric object interacting with one or several confined 
electromagnetic modes. We discuss the assumptions that 
are taken and derive the complete Hamiltonian. 



A. Assumptions 

In this description of the interaction between a dielec- 
tric with a center-of-mass position f and a photonic field, 
the following assumptions are taken: 

1. The object has a volume V, a density distribution 
p, and a mass M — pV. Note that the density dis- 
tribution is assumed to be homogeneous for sim- 
plicity. In contrast to the common assumption, see 
e.g. [53], we do not restrict the size of the dielectric 
to the sub-wavelength scale of the light field, but 
allow for arbitrary sizes. 

2. The dielectric constant e r is assumed to be homoge- 
neous. The permeability of the object p is chosen 
to be equal to the vacuum permeability, p = po, 
which is a good approximation for the dielectric 
objects we are mainly interested in. 

3. As shown in [23], the center-of-mass (cm) mode of 
dielectrics at the micron-scale is decoupled from the 
vibrational ones. Hence, we will only consider the 
motion of the cm degree of freedom f and neglect 
its coupling to vibrational modes. 

4. We assume the dielectric constant of the object to 
be real, i.e., no absorption effects are taken into 
account. In the language of scattering theory, this 
signifies that only elastic scattering processes are 
accounted for. 



Sec. |n} 

we assume the electromagnetic field to be 
scalar and neglect polarizations for a better read- 
ability of the equations. The derivations for polar- 
izations can be carried out in full analogy. We use 
the results including polarizations in the analysis 
of cavity optomechanics with levitating spheres in 
SeclVj 

6. We assume that all photons are scattered into the 
bath modes. This is a valid assumption for geome- 
tries that do not fit the cavity's geometry like, e.g., 
spheres, whereas for membranes the scattering into 
the cavity mode has to be taken into account 43J. 



B. Hamiltonian 

The Hamiltonian consists of three parts, 

H to t = Hu + Hl + i?LM : (1) 

the cm motion of the dielectric is described by Hm, the 
energy of the electromagnetic field by H^ , and the inter- 
action between the light and matter is given by -Hlm- For 
the master equation description that we want to pursue 
in the proceeding, it is useful to divide the total Hamil- 
tonian into 



H to t — H s + Hb + H 



BS ■ 



(2) 



where Hg denotes the Hamiltonian describing the system, 
Hb denotes the part describing the bath and Hbs the 
coupling between the two. Each of these terms will be 
defined in the following. 



1. The kinetic energy 

The motion of the free untrapped dielectric is described 
by Hm — P 2 /(2M), where p denotes the momentum op- 
erator of the cm coordinates in the direction we are inter- 
ested in. While the dielectric object we investigate may 
have an arbitrary three-dimensional shape, we consider 
only its motion in one dimension. Due to the harmonicity 
of the trap, the coupling between the different directions 
can be neglected. Nevertheless, in many cases it might 
still be necessary to control the motion in the other direc- 
tions, e.g., via feedback cooling [44] . In particular for lin- 
ear or quadratic potentials, also the coupling to internal 
vibrational modes of the sphere can be neglected. In the 
absence of an additional external potential, the Hamilto- 
nian can be diagonalized in the basis of the vibrational 
cigenmodes. Adding an external potential leads to some 
coupling between the cm degree of freedom and the vi- 
brational modes. The frequency of the vibrational eigen- 
modes is roughly given by w„ cx n c sound /R, where c SOU nd 
denotes the sound velocity and R the extension of the 
dielectric. For micron-scale objects it is several orders of 
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magnitude larger than the trapping frequencies typically 
achieved for the cm degree of freedom. This enables one 
to adiabatically eliminate the vibrational modes merely 
leading to a negligible renormalization of the system's 
energy. A detailed discussion using a theory of quantum 
elasticity can be found in [23] . 



2. The energy of the free electromagnetic field 



by 



The energy of the free electromagnetic field is described 



2 MX 



2 ( 

tot I 



4 2 ot(x) 

MO 



(3) 



where Ptot(x) is the object's polarization and the inte- 
gration is performed over the volume of the dielectric V 
with center-of-mass coordinate f. Assuming P tot (x) = 
a p -Etot( x ) and comparing the resulting relation between 
the polarization and the electric field for the macro- 
scopic [45] and microscopic case (see [23] for a concise 
derivation), one obtains a p ~ e$e r and 



fcfo 
2 



V(f) 



dx£ , tot (x) 2 



(7) 



where e c = 3(e r — l)/(e r + 2) is defined in terms of the 
relative dielectric constant e r . Here, the cm is treated as 
an operator, such that Eq. |7| gives the coupling terms 
between the object's position and the light field. 



where eo denotes the vacuum permittivity, E to t the elec- 
tric field and Stot the magnetic one. The total elec- 
tromagnetic field can be divided into a part containing 
the continuous modes and one or several confined modes. 
The continuous part is defined as 



£ B (x) = 



„ — ikx 



H.c), (4) 



where the label B signifies that this continuum of plane- 
wave modes will generally be treated as a bath. The dif- 
ferent modes are characterized by the annihilation (cre- 
ation) operators &k (a k ) with a mode frequency uj-^ and 
a wave vector k, where we will denote k — |k|. Note that 
we set h = 1 throughout the manuscript. In the next 
step we define a confined mode of the electromagnetic 
field with annihilation (creation) operator o,q (aj), mode 
frequency loq, mode volume Vq and a mode profile given 
by /(x). Typically it describes a mode in a cavity subject 
to some boundary conditions. We label this inhomoge- 
neous part of the electromagnetic field S (for system) , it 
is given by 



£ S (x) =i 



Up 

2e V c 



(a /(x) - H.c.) . 



(5) 



The extension to several inhomogeneous modes can be 
achieved in an analogous fashion. 



3. Light-matter interaction 

The most interesting part of the Hamiltonian describes 
the interaction between the dielectric and the electromag- 
netic field. The response of the object's polarization is 
assumed to be linear to the electric field, which is ful- 
filled for the typical light intensities considered in this 
manuscript. The interaction Hamiltonian is given by 



H 



LM 



V(f) 



dxPtot(x)-Etot(x), 



(6) 



Before describing the different contributions in detail, 
we reconsider the inhomogeneous mode E$ that has been 
separated from the continuum, see Eq. (JsJ . It describes 
one (or several) mode(s) that differs from the continuum. 
While in the specific setup of optomechanics with levi- 
tating spheres both the tweezer and the cavity field con- 
tribute, we describe this mode in general as the system 
mode. Due to the high photonic occupation numbers that 
might occur in the presence of a strong driving field, it 
can be divided into a classical part and a quantum part 
by displacing the operators So = (&q) + a' (note that we 
will omit the prime hereafter). This yields an additional 
contribution to the electromagnetic field given by 



£s(x) 



_UJ0_ 

2e a Vc 



(a/(x) - H.c.) 



(8) 



where £s (x, t) is not an operator and describes the classi- 
cal part of the light field with a — (ao), the square root of 
the photon number. Plugging E tot (x) — E s (x)+£ s (x) + 
Eb (x) into Eq. (|7|) leads to different contributions in the 
Hamiltonian H to t of Eq. (JTJ). The Hamiltonian describ- 
ing the system consisting of the inhomogeneous mode 
and the mechanical degree of freedom can be written as 



-Hs = r conakan 

2M 2 



V(f) 



dx £ s (x)+£ S (x) 



(9) 



The energy of the bath modes is given by 



(10) 



where W(f) describes the interaction between different 
bath modes induced by the presence of the dielectric, 



W(f) 



e c £o 



V(f) 



dx££(x). 



(11) 



The interaction between the system and the bath modes 
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is denoted by 



in Sec. UTTCl 



H B s = -e c eo / cbc (f s (x) + E S ( X )) E B {x). (12) 
JV(f) v ' 

The noninteracting part of the Hamiltonian describing 
the energy of the system and the bath is given by 



H o — H s + Hb ■ 



(13) 



In quantum optics, the interaction between different bath 
modes, Eq. (Ill, is commonly neglected. In contrast, 



when describing the scattering of light from larger ob- 
jects, interactions among the bath modes have to be 
taken into account, such that it is no longer justified to 
neglect W(f). H thus effects a coupling between differ- 
ent modes of the bath, such that the bath operators are 
not a diagonal basis anymore. We demonstrate in the fol- 
lowing how this problem can be addressed and connected 
to a description within scattering theory. 



III. MASTER EQUATION FOR ARBITRARY 
DIELECTRICS 



We give a concise description of the two modes of the 
system we are interested in, the mechanical mode de- 
scribing the center-of-mass motion of the dielectric and 
the cavity mode of the light. Therefore, we trace out the 
other modes of the electromagnetic field, the free modes. 
The typical quantum-optical approach to these systems 
is the method of Born-Markov master equations, where 
the bath is eliminated to derive a description exclusively 
for the system's dynamics. The Hamiltonian is split into 
a part describing the energy of the system and the bath, 
Ho, and the interaction between the two, Hbs- For typi- 
cal quantum-optical systems, Ho is diagonal in the bath 
operators Sk as interactions among them are negligible, 
such that the transformation to the interaction picture is 
straightforward. The difficulty we confront when describ- 
ing the interaction between light and a dielectric sphere 
in a cavity larger than the wavelength is that due to the 
large number of scattered photons, interactions within 
the bath, given by Eq. (11), have to be taken into ac- 
count. This effects a Hamiltonian which is non-diagonal 
in the bath operators Sk- 

The strategy to approach this problem is to first solve 
the equations of motion, effected by the interaction with 
the dielectric, for the bath operators. Connecting these 
expressions to the Lippmann-Schwinger equation for the 
scattering process of a single photon, we give the solu- 
tion for the bath operators containing the full scattering 
interaction in Sec. |HIA} S ubsequently, we derive the mas- 
ter equation in Sec. Ill B| describing the cavity mode and 
the center-of-mass mode in this new basis of bath oper- 
ators, enabling one to express all quantities in terms of 
scattering operators. Finally, we specify this approach 
to optomechanical systems, assuming a strongly-driven 
cavity and the Lamb-Dicke regime for the cm operator 



Solution of the scattering equations for the free 
field 



The equations of motion of the electromagnetic field 
in the Heisenberg picture are determined and connected 
to the Lippmann-Schwinger equation. We are only in- 
terested in the homogeneous part of the electromagnetic 
field Eb(x, t) given by Eq. Q and assume that no in- 
homogeneity (i.e., cavity) is present, leaving us with a 
system fully described by H B , Eq. (10). We keep the 



center-of-mass operator r in the equations of motion, but 
neglect its action for now assuming M ~ > oo. The Heisen- 
berg eqs. of motion can thus be determined as 



d k (t) = -iu k Qk(t)- Ul Zn' ^ / dxE B {x,t)e ikx . 
V 2(2vr) J J v{f) 



(14) 



Let us first define E B = E ( B +) +E B \ where E ( B } = E ( B +)t 
with 



(2tt) 3 / 2 

and the incoming field is given by 



(/k(/ ^k e -ikx- kW (15) 

ze 



4"£(m) 



./k < /^e- ikx a k (0)e- iWk *. (16) 
2e 



(2tt) 3 / 2 

With these definitions at hand, we can close the set of 



equations given by Eq. ( 14 ) by carrying out the following 



steps: we formally integrate Eq. (14) over time, multiply 



both sides by i-\/wk/2eo(27r) 3 e lkx , and take the integra- 
tion over k to obtain 

4 +) (x,f,t) =E$l(x,t) 

v(f) J 2(2^)3 



+ E^~ ) (x',t) / dre i (" 0+Wfc)T 

./o 



(17) 



Note that we assume the spectral distribution of the elec- 
tromagnetic field to be peaked at a certain frequency 
wo and we thus have introduced the slowly-varying field 
E' B (x, t) here, ^(x, t) = e ±iuot ^(x). This justifies the 
assumption that E'jf (x, t) remains constant on the time 
scales of the system's evolution which allows one to take 
it out of the integration in Eq. (17). Integrating dk in 



Eq. (17) yields a function that decays quickly in r. This 



allows for an extension of the upper integration boundary 
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t to oo (Markov approximation) and hence yields 

dx' / dk 

K^r Jv(f) J 
iklx . , — ^ — 4+v,i), 



w k - w + 17 



(18) 



where the limit 7 — > + is understood. Taking the inverse 
transformation, the field operators can be written as 



Sk(0 =a k (0)e 



— lUJ^t 



v(t) 



dx' / dk'x 



(19) 



W k ' - wo + 17 



This equation for the operators of the electromagnetic 
field resembles the Lippmann-Schwinger equation [461 
Wf\ . In order to connect the two descriptions, enabling 
one to employ solutions known from scattering theory in 
our approach, we proceed in the same way defining 

V k ,k'(r) = ^/ dx'VZJk^ke^'^' (20) 

1 JV(f) 

as the matrix elements of the operator describing the 
scattering interaction. In analogy we define the transition 
matrix 7k,k'(^), given by 



T k ,k'(r) = J dk"V k ,k»(r)x 



x <S(k' - k") + 



Tk» k >(r) 



(21) 



Note that both Vk,k' (?) and 7k, k' (?) are operators for the 
center-of-mass degree of freedom but not for the photonic 
ones. That is, if we fix f, neglecting the object's motion, 
V^k'^) and Tk,k>(r) are simply numbers without any 
operator-character. By iteration, the transition matrix 
describes scattering processes to all orders of perturba- 
tion theory, as illustrated in Fig. [T] (an explicit formula 
for the expansion of 7k,k'(^) in terms of Vk,k'(^) can be 
found in standard textbooks [4 7) ) . 



Subsequently, Eq. (21) enables one to rewrite the total 



time evolution of the operators, Eq. (191, as 



a k (t) = / dk'e 



<5(k-k') 



T k ,k'(r) 



(22) 



Ok'(0). 



This expression is equivalent to the classical field equa- 
tions given by Eq. (Al). For its solution we can thus 



rely on the variety of methods that has been developed 
during the past decades described in Sec. |XJ 

A useful relation that will be used to simplify the com- 
putation of transition amplitudes for spheres in Sec. [VJ 
is the optical theorem connecting the scattering ampli- 



k,k n (r = 



w -. .. 

/_ x 



✓ 'X U* X 



FIG. 1: Graphical illustration of the transition matrix 
7k,k if) as an infinite series in the light-matter interaction. 
The first term denotes the direct interaction between the di- 
electric and light, the second one a process, where a photon 
is virtually absorbed and reemitted, the third one a process, 
where two intermediate photons are involved, etc. 



tude in forward-direction to the scattering in all other 
directions 146 : 



9[7k,k'(r) 



dk'|Tk,k<(r)| 2 <5K-WkO- (23) 



Before we continue the analysis, to ease the notation, 
it is useful to define the space of mode functions in which 
the matrices T~\ty(r) and Vk,k'(^) act, and consider them 
as operators, i.e., Vk,k'(^) = (k|V r (?')|k'), where |k) are 
the basis vectors of such a space. They can be viewed 
as mode functions with momentum k. As we describe 
scattering out of the cavity mode in this article, let us 
now define the transition amplitudes for mode shapes 
different from plane waves and express them in the basis 
|k). For the Born approximation of scattering theory [62], 
consisting in setting T(f) ss V(r), we obtain 



V k ,c(r)= / dk'(k|f(f)|k')(k'|c) 



(24) 



where |c) describes the mode function of the cavity, 
which can be written as (x|c) = f(x)/y/Vc in position- 
representation, where /(x) is assumed to be real. Evalu- 
ating this expression yields 



Vk,c(r) 



4V C 



dx/(x)e i 



kx 



(25) 



V(r) 



where we have used that the distribution of is peaked 
around ujq, allowing for the substitution Wk ~ wo- 
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B. General master equation for the cavity and the 
center-of-mass mode 

Within the Born-Markov approximation, the master 
equation describing the system's full dynamics is given 
by 



Ps(t) 



/ tr B [i% s (t),[HiJ S (f 
Jo 



r),ps(t) ® p B ]]dr. 

(26) 

The Born-Markov approximation consists in the follow- 
ing assumptions: the density matrices of the system 
and the environment are considered to be separable, 
Ptot = Ps ® Pbi and correlations between bath opera- 
tors are taken to decay quickly. Furthermore, the bath 
is assumed to remain unchanged during the interaction 
with the system, /5b (t) ps Pb(0). This is valid given that 
the bath is very large and the effect of the interaction 
with the system can be neglected. 

Moreover, for typical quantum-optical systems interac- 
tions between different bath operators a k ,a k , , Eq. (Ill, 
are negligible, i.e., W(f) ps 0, such that 



H pa H s + I dkuj k ala k 



(27) 



While these approximations are typically fulfilled for 
point-particles, difficulties are encountered when extend- 
ing the method to larger objects. It is in particular the 
negligence of interactions between different bath opera- 
tors that is no longer justified. More specifically, we re- 
alize that contributions, where interactions among bath 
operators are taken into account to different orders, scale 
as oc (R/X) 2n . Here, R denotes the dimensions of the ob- 
ject, A the wavelength of the inhomogeneous light mode, 
and n the nth order of the multiple scattering process. 
Consequently, an approach where these interactions are 
accounted for is necessary. This is done by including 
the correlations between the bath operators described by 
W(f) in the Hamiltonian that is used to transform to the 
interaction picture, 



#bs(*) = e iHot H BS (0)e 



-Wat 



(28) 



To find a solution, we connect this approach to the de- 
scription within scattering theory given in Sec. |III A| 
Based on this analysis, we can develop a master equa- 
tion that accounts for interactions among different bath 
modes. As an example we now discuss the first term of 
the master equation, where all operators are in front of 



the density matrix, in more detail: 

<h I dke iw ° T (O|JX*>r)a k (0)|O)x 



Ps 



%a 



o j (29) 
x (n\a k (0)T(t-T,f)\n)p s + ... 
Here, |0) denotes the vacuum state and we have defined 



dx/(x)£ B (x,i), 



(30) 



V(f) 



where counter-rotating terms have been neglected. Let 
us now connect Eq. (29) to the description in terms of 



mode functions in the scattering picture. First, we shift 
the time dependance of F(t, f) to the operators by 



(n\F(t,r)at(0)\n) = (Sl\T(0,f)a{(-t)\n), 



(31) 



where the invariance of the vacuum state under time 
evolution has been used. In order to make the proce- 
dure more transparent, as a first step, only the 0th or- 
der Born approximation of scattering theory is identified. 
Subsequently, the treatment is extended to a description 
of all orders. The lowest order of the Born series gives 
&k(i) ~ Ok(0) cxp(— icJkt), so that we need to evaluate 



(01^(0,^(0)10) = ^^ 



V(r) 



dxf(x)e 



— ikx 



(32) 



Recalling the definition of the expectation value Vk,c(^) 
in the scattering picture, Eq. (25), we identify 



(o|j-(o,f)4(o)|o) = v; k (f). 



(33) 



The same procedure can now be applied without taking 
the Born approximation and considering the full transi- 
tion matrix, by plugging Eq. (22) into Eq. (31), which 
yields 



(O|J-(0,r)4(-« 

= / dk'V c *. k ,(f) (V-k' 



T k *, k (f) 

W k - 0JQ - 17 



(34) 



where Eq. (21) has been used. 



All other terms of the master equation can be deter- 
mined in full analogy yielding 



Ps =i\ftS,H s + Hm] + 



J dkS(u 



w k - w ) I 27 k , c (f : )ao/5sa '^tk(^) ~ |7k,c( 



| 2 a a ,/5s 



(35) 
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where Hg is the system Hamiltonian given by Eq. {§) 
and H r „ the renormalization 



dk\% 



V - 



(36) 



where V denotes Cauchy's principal value. Note that 
a similar master equation for the cm degree of freedom 
has been discussed in the context of scattering of air 
molecules 



C. Master equation for the optomechanical setup 



In this section, we adapt the general master equation 



Eq. ( 35 ) to the specific optomechanical setup we are in- 
terested in. Therefore, we take the following approxima- 
tions: 

1. The inhomogeneous mode is assumed to be a 
strongly-driven cavity effecting large cavity occu- 



transforms and the analytic property of the function 
7^ c (r)|f=o7^* k (0), see App.[C]for an explicit analysis. 

The only contributing terms are oc f 2 and describe a 
renormalization of the trapping frequency of the dielec- 
tric provided by the optical tweezers and decoherence of 
the cm operator. The renormalization of the trapping 
frequency ui t = u t + A M can be simplified exploiting the 
analytic properties of the functions and is given by 

A M = \a\ 2 xl J dVP^^ [T^ c (f)T c l(f)]. =0 . (37) 

The decoherence of the mechanical motion is described 
by 



£ M [p s ] = r ( 2(8 + P)p s (b + St) - [(6 + 5t) a ,^g 



with 



strongiy-arivcn cavity cnccmig large cavity occu- r 

pation numbers n phot = |a| 2 , such that a > (a ). r = \a\ 2 %l J dkS{u u - uj ) [K iC (r)T£(r)] f=Q 



(38) 



(39) 



This enables one to neglect certain terms in the 
master equation. 

2. We assume the Lamb-Dicke regime: the dielectric is 
positioned close to the maximal slope of the stand- 
ing wave in the cavity and close to the minimum 
of the harmonic trapping potential of the optical 
tweezers. The motion around its cm position is con- 
sidered to be small, such that the Lamb-Dicke pa- 
rameter rj = kAf <C 1 (with Af = yj (f 2 ) — (r) 2 ), 
facilitating an expansion of the transition operator 
matrix elements 7k, c (^) m k^- 

Displacing the cavity operator by a such that ~~ ► «o+oj 
and expanding the transition operator to second order, 
7L,c(0 ~ 7k, c (0) + U c {r)\?= r + T£ >c {r)\ r = Q r 2 leads to 
a master equation, where we take into account terms 
that are at most of quadratic order in the cavity opera- 
tors fio,aJ and the cm operators r = xo(b + Here, 
Tk c (?) = d n %L,c{r) / dr n denotes the nth partial deriva- 
tive and Xq = \flj2MuJt the zero-point motion of the 
center-of mass mode. 

In the following we give an interpretation of the dif- 
ferent contributions to the master equation and indicate 
which terms yield a renormalization to the Hamiltonian, 
can be neglected, or describe decoherence. We describe 
these terms in decreasing order in a. 



Contributions oc \a\ 



The decoherence of the cm thus depends on the form 
of the transition amplitudes with respect to the cm posi- 
tion. The physical process underlying this effect is recoil 
heating via photon scattering. 



2. Contributions oc a 

Also for contributions oc a*ao|7k, c (0)| 2 , the analyticity 
of the transition operator can be exploited. Applying a 
Hilbcrt transformation, we can show that these contribu- 
tions cancel, see App.[C]for a more detailed analysis. 

Terms oc a*dof effect both a coherent and an inco- 
herent contribution. The incoherent part describes de- 
coherence of the mechanical and the light degree of free- 
dom and can be shown to be negligible as demonstrated 
in App. [C] In contrast, the coherent contribution yields 
a non-negligible renormalization of the optomechanical 
coupling g — g + g rn defined by 



9r 



a xq 



dkV- 



T k ,c(0)T c '; k (f)|, =o . (40) 



Furthermore, terms oc f 2 describe decoherence of both 
the mechanical mode and the light mode. Comparing 
to Eq. (39 1 for the cm mode, these contributions are 



suppressed by 1/a and can thus be neglected. Also for 
the cavity mode, these terms are negligible, given that 
rr 2 < 1/a. 



The largest contribution to the master equation are 
terms oc |a| 2 |7k, c (0)| 2 . As they do not contain operator- 
character they cancel. 

The next order in the Lamb-Dicke parameter 77 is given 
by terms oc f, which can be shown to vanish using Hilbert 



3. Contributions oc a ao 

Terms oc iTk^ (0) | 2 yield both a coherent and an inco- 
herent contribution describing a renormalization of the 



resonance frequency of the cavity and a part describing 
the cavity's decay: 



'[ps] =« (2a a p s al - [aja , ps]+"j (41) 



with 



k = J dk<5K- Wo )|T k , c (0)| 2 



(42) 



The renormalization of the cavity's resonance frequency 
ljq = ljq + A L is defined by 



A L = / dkV 



1 



■|Tc, k (0)| 2 



(43) 



These are the only non-vanishing contributions as terms 
oc f are suppressed by the Lamb-Dicke parameter r\ and 
terms oc f 2 even by i] 2 compared to Eqs. (41), (43). 



4- Final master equation 

To summarize, we identify the contributions to the fi- 
nal master equation: 



(44) 



Ps =i[ps, H S + H m ] + C M [p s ] + C L [p s ] 
They can be grouped as follows: 

1. Contributions of Hamiltonian-type, 

H s + H m = - Sala + Cb t b r b + g(a a + aJ)(S + 6 f ), (45) 

where the frequencies and couplings stemming from 
the system's Hamiltonian i?g, given by Eq. @, are 
renormalized by 

H xn = A M tfi> + A L d a + g ln (a + a$)(S + S+). (46) 

The corresponding renormalizations arc defined by 
Eqs. j37|, d40b, pSb. Note that the Hamiltonian of 



about the scattered fields can be obtained by applying 
the quantum regression theorem. The scattered light is 
directly accessible in experiments and can, e.g., be used 
to monitor the cooling of the mechanical motion [53] - 
To complement the analysis given here, we show how to 
derive the scattered fields directly in App. [E] within an 
approach similar to Wigner-Weisskopf, but accounting 
for interaction processes between the bath modes. 



IV. COOLING 

Before applying the master equation discussed in the 
previous section to the particular case of cavity optome- 
chanics with optically levitating spheres, we provide a 
general description of optomechanical cooling and the 
minimal phonon number attainable with master equa- 
tions of at most quadratic order in the operators of the 
mechanical and the cavity mode. Cooling is in gen- 
eral a vital ingredient in any attempt to demonstrate 
quantum-mechanical behavior. Most descriptions make 
certain approximations to ease calculations, e.g., the 
sideband regime is commonly employed in optomechani- 
cal setups [301 OS] • Here, the system will be treated in the 
most general way not relying on any approximations (as 
some of them might not be fulfilled for larger objects). 
The master equations we are interested in, are typically 
of the form given by Eq. ( 44 ) . 



The mean phonon number (wb) is coupled to all 
other expectation values of combinations of the operators 
do, a ,b, and yielding the Eqs. of motion in matrix 
form 

v = Mv + c, (47) 

where 

v= ((tib),(a tf),(aiai),(aib),....) T , (48) 
M denotes the interaction matrix, and c a constant vec- 



Eq. (45) has been transformed to a frame rotating tor. The master equation Eq. (|44j) keeps the Gaussian 



at the laser frequency u>l, where S now denotes its 
detuning from the cavity resonance frequency ojq. 

2. The recoil heating via photon scattering of the cm 
mode yields £ M [ps], given by Eq. (38). 



character for an initially Gaussian state. Subsequently, 



also the system of equations Eq. ( 47 1 is Gaussian in the 



operators ao,al... and linear in the expectation values 



(w b). Since Eq. (47) represents a closed system of equa- 
tions, they can be solved as 



3. The decay of the cavity mode due to the presence of 
the object yielding £ L [ps], is described by Eq. (41 ). 



Consequently, all frequencies, couplings, and decay rates 
are renormalized taking into account all terms beyond 
the first Born approximation of scattering theory. This 
enables one to use exact solutions if available, or in gen- 
eral to truncate the perturbation series in a controlled 
way. 

While this master equation only contains the time evo- 
lution of the cavity and the cm operators, information 



yielding 



Mo) 



v = —M~ 



(49) 



(50) 



for the steady state. The steady state phonon number 
can be extracted from this quantity as 



n = lim Woo (b t b) 



Ai - A 2 



A, 



A a 



(51) 
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with A 1 = 32g A S[AS 2 K + K, 3 + 16S(T-K,)Cb t + 8Kut], A 2 = 
-T[4S 2 +n 2 ]Cb t [16S 4 +8S 2 (k 2 - 4£ 2 ) + (k 2 + 4w t 2 ) 2 ], A 3 = 
4g 2 {-K(4<5 2 + K 2 ) [k 2 + 4(6 - u t ) 2 ]Cj t + 2r<5[32w t 4 + (4<5 2 + 
k 2 ) 2 + 4(-205 2 + 3k 2 )w 2 ]}, and A 4 = 64g 2 6nuj t [mg 2 6 - 
(4<5 2 + k 2 ) a) t ]- This solution might contain unphysical 
results. To verify that n is indeed a steady state of the 
system, the eigenvalues of M additionally have to fulfill 
3?[eig[M]] < 0. In general, all parameters are determined 
by the properties of the system, solely the detuning 5 
can be chosen. According to the definition of Eq. (451, 
8 > denotes red detuning and 6 < blue detuning. To 
obtain the optimal point for cooling, one consequently 
has to optimize n with respect to 8. Let us now compare 
this exact solution to the one obtained after an adiabatic 
elimination of the cavity mode. Starting from Eq. ( 44 1 , 



we eliminate the cavity mode assuming that its decay 
rate is much larger than the coupling between the me- 
chanical degree of freedom and the light, n 3> <?. In 
this case it is justified to assume that the cavity is ei- 
ther empty or contains only one photon, therefore reduc- 
ing the master equation to the one-excitation manifold, 
described by poo, Pw, Poi, Pn- Due to the fast decay of 
the cavity mode described by k, the change of all con- 
tributions involving an excitation is approximately zero, 
finally yielding an equation of motion for the empty cav- 
ity poo- After carrying out a rotating wave approximation 
assuming Co t ^> \g 2 / (k + i(8 ± £j t ))\ , the final steady state 
phonon occupation is given by: 



[4g 2 K + T(k 2 + 4(5 + £ t ) 2 )] [n 2 + 4( 5 - d> t ) 21 



^-adiab 



To obtain the minimal occupation number, this equation 
needs to be optimized with respect to S. Comparing the 
adiabatically-eliminated solution to the exact one, it be- 
comes clear that the approximation breaks down in the 
strong-coupling regime g m uj t , where the rotating- wave 
approximation is no longer valid, see Fig. [2] for an illus- 
tration. Note that Eq. (52 1 can be derived from Eq. (51) 



taking the approximation k» j. In case we choose the 
detuning 6 = Cjt and u>t 3> n, we obtain the minimal 
occupation number in the sideband regime 



4£j t ) 



1 

4C : 



(53) 



where C denotes the cooperativity, given by Eq. ( 62 1 . 



CAVITY QUANTUM-OPTOMECHANICS 
WITH LEVITATING SPHERES 



The theory developed in the previous sections provides 
a general framework to describe dielectrics interacting 
with one specific mode of the electromagnetic field. Here 
it will be used to describe cavity-optomechanics with 
levitating spheres [221 [2H H5J- Due to the levitation, 
their mechanical degree of freedom (i.e., the motion of 




0.5 1.0 1.5 2.0 %.0 0.5 1.0 1.5 2.0 

6/u>t S/u>t 
80 r. 





FIG. 2: Minimal phonon numbers for different detunings com- 
paring the exact result (blue straight line) to the one with 
the cavity mode eliminated (red dashed line). Differences be- 
tween the two solutions evolve as g is increased. Only the 
regions, where a steady state is attainable have been plot- 
ted, they all lie within the red-sideband regime. For blue 
detuning or an optomechanical coupling which is too strong, 
the system is heated and no steady state can be obtained. 
Upper left pannel: Sideband-resolved regime with weak cou- 
pling, k = 0.3 Wt,r = 0.03 ait,£f = 0.07 ti t , Upper right pan- 
nel: k = 0.3 cDt,r = 0.03 &t,g ~ 0.3 ait, Lower left pan- 
nel: Bad-cavity limit for weak coupling k = 3 iit,r = 0, g = 
0.1 il>t, Lower right pannel: Bad-cavity limit for strong cou- 
pling, K = 3iit,r = 0,j = 0.864 cD t 



their cm) is prevented from clamping losses, which are 
the main source of decoherence in most optomechani- 
cal systems. Hence, they hold the promise for a variety 
of applications, ranging from the efficient implementa- 
tion of protocols to realize non-Gaussian states [23] to 
the preparation of large superpositions of their cm po- 
sition [T?l [T5]. In particular this enticing perspective 
has the potential to realize tests of quantum mechan- 
ics in an entirely new parameter regime. All of these 
results are only valid for nano-objects smaller than the 
optical wavelength, enabling one to neglect multiple- 
scattering effects. Although schemes to Doppler-cool 
dielectric spheres using Mie resonances have been dis- 
cussed |54) . cavity-optomechanics with larger levitating 
objects has not been described before. In particular in 
the light of recent experiments on feedback-cooling of 
a microsphere [36 , a description of this new parameter 
regime is timely 

We thus proceed to determine the optomechanical pa- 
rameters and final occupation numbers for levitating di- 
electric spheres. The setup is similar to the one discussed 
in [53], but is extended to a description of spheres larger 
than the cavity wavelength. It is sketched in Fig. [3] 
a classical light field, effected by a retro-reflected opti- 
cal tweezer, £ tw (x, t), creates a harmonic trap for the 
cm of the dielectric, (note that trapping via a strongly- 
populated cavity mode can be described in full analogy) . 
Besides, a second cavity field E cav (x, t) is used to manip- 
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FIG. 3: Schematic representation of the setup: A dielectric 
sphere of a radius similar or larger than the cavity wave length 
is trapped by optical tweezers providing a trapping frequency 
ujt. It is placed inside an optical cavity, where a second laser 
is used to optically manipulate the dielectric's center-of-mass 
degree of freedom. 



ulatc it, such that 



#sCM) = E cav (x,t) 

£sCM) = £tw(x,f) +£ca.v(x,t). 



(54) 



The optical tweezer used for the trapping is given by 



£tw(x,t) = E 



W t 
W(y) 



exp 



W(yf 



(55) 



where E = [P t j ' (eoCKW 2 )} 1 ! 2 , P t is the laser power, W t 
is the laser beam waist, W[y) = W t [l + (y\/ '(irW?)) 2 } 1 / 2 
and we assume the beam to be aligned as sketched in 
Fig. [3j While we are only interested in the classical part 
of this field as it is used solely for the trapping, we include 
both the quantum and the classical part of the cavity 
field consisting of a standing wave in z-direction and a 
Gaussian profile in x- and y-direction, 



w 



E cav (x,t) = L — fz (/cav(x)a - H.c.) 



w 



(56) 



£ cgcv (x,t) = i,/— — (/ ca v(x)a - H.c). 



This equation denotes the cavity field in the displaced 
form, where \a\ 2 — n p h is the mean number of photons 
in the steady state, n p h = y /, 2P c n/uJo/(i5 + k), with P c 
being the power of the driving laser. The mode function 

is given by / Gav (x) = exp ( — ^^gf ) cos(k 0tZ z-ip), where 

ip denotes the equilibrium position of the dielectric, ko 
the wave vector of the cavity light, and the mode volume 
is given by Vq — LttWq /4 with L being the cavity length 
and Wq its waist. While the classical term merely yields 
a shift of the trapping frequency and the equilibrium po- 
sition, the quantum part of the mode function is used 
to manipulate the cm degree of freedom of the dielectric 
including the part describing the opto-mechanical cou- 



pling. Note that we only consider one mode of the cavity 
here, higher harmonics supported by the cavity are not 
included, they are contained in the continuum of homo- 
geneous modes and coupling to them is treated as losses. 
In case of using a second cavity mode for the trapping 
instead of the tweezer, Eq. ( 56 ) simply has to be summed 



over several modes with different profiles. 

The full dynamics of the system is obtained taking into 
account the coupling of the tweezers and the cavity mode 
to the vacuum modes, given by Eq. (4|. The full master 
equation is described by Eq. ( 44 ) with the corresponding 



decay rates given by Eq. ( 39 ) and Eq. ( 42 ) , where T con 



tains contributions of the cavity mode and the tweezers. 

The specific description of spheres is eased by the avail- 
ability of an analytical solution, the Mie solution, for the 
scattered fields and cross sections [HJ 22] • The Mie solu- 
tion is based on expanding the incoming electromagnetic 
field in spherical waves. This expansion suits the sphere's 
geometry and it is thus possible to apply boundary con- 
ditions to determine the scattered fields. Note that while 
the polarization of the electromagnetic field has been ne- 
glected to ease the notation in the previous sections, we 
take it into account here. The Mie solution is defined for 
plane-wave states, and we use the relation 



7k, k' = 



27ra;k 



/(k,k') 



(57) 



to formulate the solution for the cavity field in terms of 
classical amplitudes. 

For spherical objects, assuming a vanishing absorption, 
$s[e r ] s» 0, it is possible to connect all quantities to the 
classical scattering amplitude in the forward-direction us- 
ing the optical theorem, Eq. ( 23 ) , which yields 



/(k,k) = ^^(2n + l)(a n + 6 n ). 



(58) 



n=l 



The coefficients a n ,b n depend on the dielectric constant 
and the radius of the sphere, they are defined in terms of 
spherical Bessel functions. We refer the reader to stan- 
dard textbooks [HJ 02] for their specific form. Hence, 
the optomechanical parameters can be determined: 

1. The optomechanical coupling is defined by 

XqTTC 



9 =ot 



2k V c 



sin(2<p)x 



x ^ 



(59) 



where c denotes the velocity of light and tp the po- 
sition of the sphere in the cavity. 



2. The total cavity decay rate is defined by K tot = 
kq + k, where kq is the intrinsic cavity decay rate 
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resulting from imperfections in the mirrors and 



2klV c 



3?[^(2n + l)x 



(60) 



x (l + (-l)"cos(2<p))(a„ + &„)]. 



3. The recoil heating of the dielectric due to scattering 
of cavity photons can be computed as 



jicav 



K[> (2n + l)x 



V c ^ ' (61) 
x (l + (-l)"cos(2^))(a n + &„)]. 



Note that the recoil heating from the trapping 
lasers can be obtained in full analogy by inserting 
the tweezer mode. 

Besides the minimal phonon number n m - ln describing the 
possibility to cool the system (close) to its quantum- 
mechanical ground state, another figure of merit to de- 
scribe the cavity-optomechanical properties is the coop- 
erativity C. This measure for the coherent coupling be- 
tween the motion and light is defined by 



Tk 



(62) 



and depends on the size of the particle and its position 
in the cavity. It is in particular essential to have a suffi- 
ciently high cooperativity to perform protocols coupling 
the cm to the light [23]. Assuming that the object is po- 
sitioned at the maximal slope of the standing wave and is 
much smaller than the laser's waist, the asymptotic form 
of the cooperativity is given by 



C oc 



l/e 2 c k%R 6 ifi?< A; 
l/k 2 R 2 ifi?>A 



(63) 



under the assumption that the laser's waist is larger than 
the object. In case the dielectric is not fully covered by 
the laser's waist, the beam's Gaussian shape has to be 
taken into account {55) leading to an even lower coop- 
erativity. In the proceeding, both the minimal phonon 
number and the cooperativity are used to quantify the 
system's performance as an optomechanical setup. The 
optomechanical parameters are thus determined as illus- 
trated in Fig. [4] for varying sphere sizes. The experimen- 
tal parameters are chosen as follows: 

• Dielectric object: We assume spheres of fused 
silica with density p — 2201 kg/m 3 , a dielectric 
constant 3ft[e r ] = 2.1 and 3[e r ] - 2.5 x 10~ 10 . We 
vary their radii between R = lOnm — 2/im and posi- 
tion them at the maximal slope of the cavity field, 
if = 7r/4. 

• Cavity: We assume a confocal high finesse cavity 
of length L — 4 mm and finesse T — 5 x 10 5 leading 
to a cavity decay rate /to = cir /2FL = 2ir x 44kHz. 




FIG. 4: Cavity-Optomechanical parameters for different sizes 
of the object: Upper left panel: Optomechnical coupling g, 
Upper right pannel: Recoil heating of the center-of-mass F, 
Lower left pannel: Cavity decay rate k, Lower right pannel: 
Cooperativity C 



This cavity is impinged by a laser of power P c 
1 mW, wavelength A = 1064 nm, which gives a 
waist of Wq = \J \d/2-K 26 (im 

• Optomechanical parameters: The tweezers are 
constructed with a laser of wavelength A = 1064 nm 
and a lense of high numerical aperture. They sup- 
ply a harmonic trap for the object of frequency 
uj t = 27r x 136 kHz in the transversal direction 
and a slightly smaller one in the direction of light 
propagation. The cavity photons have a frequency 
cj = 27r x 2.8 x 10 14 Hz and the steady state photon 
occupation is \a\ 2 « 3.7 x 10 8 . 

Considering the graphical illustration of the optome- 
chanical parameters in Fig. [3J the absolute value of the 
optomechanical coupling g first increases with the radius 
R reaching a local maximum at R sa 260 nm, then de- 
creases and even vanishes at R as 370 nm. In the fol- 
lowing it continues these oscillations. The decoherence 
rate of the cm motion first increases oc R 3 , then be- 
gins to fall off for R > 600 nm. This is due to its de- 
pendence on the ground-state size and the cross section, 
where the scattering is described by the Rayleigh cross 
section cx R 6 for small objects, to give way to a scal- 
ing oc R 2 in the limit of geometrical scattering and the 
squared ground state size, which is oc R~ 3 . Also the cav- 
ity decay rate increases oc R 6 at first, then shows some 
plateaus to finally converge to a scaling oc R 2 . Conse- 
quently, the cooperativity first decreased immensely to 
exhibit oscillations later on. These oscillations can only 
be predicted taking multiple-scattering processes into ac- 
count. Nevertheless, the maximal values of the coopera- 
tivity are merely C « 0.05. The minimal phonon number 
is obtained by minimizing the function n described in 
Eq. (51) with respect to S. While ground-state cooling 



is feasible for spheres R < 250 nm, only relatively large 
final phonon numbers can be achieved for larger spheres, 
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FIG. 5: Minimal phonon number attainable for different 

sphere sizes with the experimental parameters given in the 
main text. 



e.g., 7i m in ~ 350 for R s» 1.3/im. 



VI. CONCLUSION AND OUTLOOK 

We have derived a full quantum theory to describe 
the coupling of light to the center-of-mass motion of 
non-absorbing dielectrics of arbitrary sizes and shapes 
in optical cavities. The common approach to the de- 
scription of sub-wavelength dielectrics and light is within 
the framework of Born-Markov master equations [51)] . 
For quantum-optical systems, interactions among differ- 
ent bath modes are commonly not present or negligible. 
However, as the size of the dielectric is increased to the 
order of the cavity wavelength, this approximation is no 
longer justified. In this article, an approach to take these 
interactions into account has been developed. This is 
achieved by describing the bath operators of the electro- 
magnetic field in terms of their scattering solution in the 
presence of a dielectric and developing a master equation 
in this new basis. Based on this method, it is possible 
either to solve the problem exactly, if analytical solutions 
are available, or to truncate the series of multi-scattering 
processes in a controlled way. The resulting master equa- 
tion of Lindblad form gives a full description of the sys- 
tem's dynamics, more specifically, decoherence rates for 
both the cavity and the light mode as well as renormal- 
izations to the Hamiltonian can be calculated. 

In the second part of the paper, we apply these meth- 
ods to optically levitating spheres of arbitrary size, de- 
termining all optomechanical parameters. In particular, 
the solutions are expanded in terms of classical scatter- 
ing amplitudes for plane- wave states, demonstrating that 
the knowledge of these amplitudes is sufficient to fully 
analyze the system. In the spherical case, an analyti- 
cal solution for the scattering amplitudes, the Mie solu- 
tion, simplifies computations. Summarizing the results, 
we find that on the one side, small spheres can be cooled 
to their quantum-mechanical ground state and can be ef- 
ficiently addressed by light due to a high cooperativity. 



On the other side, as expected, ground-state cooling is 
out of reach for larger spheres and only relatively high 
phonon occupations n min sw 350 can be obtained. 

Albeit the focus of this article is on the general deriva- 
tion of the theory and the analysis has been restricted to 
transparent spheres without internal structure, we would 
like to give an outlook on two extensions of this work that 
might significantly decrease the minimal phonon num- 
bers. First, objects which are better-suited, or even tai- 
lored to the cavity mode, might perform more efficiently 
as an optomechanical device. Given the fact that they 
scatter most photons back into the cavity mode, losses 
are reduced yielding more benign conditions for optomc- 
chanics and thus lower phonon numbers and a higher 
cooperativity. Another possibility could be to develop 
new cooling schemes based on the additional use of in- 
ternal degrees of freedom of the sphere. As the sphere 
grows larger, the coupling between the cm and the rela- 
tive motion is increased and additionally displays some 
resonances. These resonances hold the potential to be 
used to cool the cm motion via coupling to the relative 
degrees of freedom, already being sparsely occupied due 
to their high frequencies. 
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Appendix A: The classical approach 

This Appendix sketches the solution of the equations 
of motion of the classical electrodynamic field, giving an 
overview of the possibilities and limitations of the de- 
scription in classical scattering theory. The approach in 
the classical case is to solve Maxwell's equations [45j . 
Neglecting polarizations yields 

B(x, t) = E in (x, t)+e c J dx'0(x', x)£(x', t), (Al) 

where E- lxl {x.,t) denotes the incoming electromagnetic 
field and £?(x',x) = |k | 2 exp(i|k ||x - x'|)/|x - x'| the 
propagator (krj being the wave vector of the incoming 
field). This self-consistent equation has been intensely 
studied in classical scattering theory and is in general 
only solvable approximately. There exist only few ge- 
ometries, like, e.g., a cylinder, a sphere, or an ellipsoid, 
where analytical solutions are tractable. In the special 
case of a spherical object, the scattered electric field can 
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be determined exactly by expanding the field in spher- 
ical waves and subsequently applying boundary condi- 
tions, yielding the so-called Mie solution [4"0Tl42] . Per- 
turbative approaches [57J ISE], based on the analytical 
solution and an extension of the treatment via distorting 
the surfaces at different points, only allow for calcula- 
tions of small perturbations. Numerical approaches like 
the discrete dipole ansatz [38] or the T-matrix method 
(see, e.g., [35] for an expository article) are applicable 
to a larger class of objects and are widely used today. 
Indeed, these approaches coincide with the analytical so- 
lution for perfectly spherical objects [35]. In the limit 
of very large objects, R 3> A, applying a ray-optics ap- 
proach immensely simplifies the calculation of forces on 
the dielectric [5j5]. A further analysis of the classical so- 
lution is beyond the scope of this article and we refer 
the reader to the literature, for example, [371 for a more 
detailed discussion. 



Once the electromagnetic field including the scattering 
is obtained, the classical radiation force is determined 
via the momentum conservation law: the force acting 
on the dielectric is the change in momentum of the EM 
field and can be determined from Maxwell's stress tensor. 
The total force on an object interacting with the EM field 
consists of the change of mechanical momentum and field 
momentum F to t = dP mcc ^/dt + dP^i^/dt with 



mcch 



dt 

d d_ 
dt ; fiold ~ dt 



= / dx(p e E + J x B) 
Jv 

/ ehc(E x B), 
Jv 



(A2) 



where B denotes the magnetic field, p e the charge den- 
sity and J the current. Rewriting and manipulating 
this equation (see [45] for details) yields a formulation 
in terms of Maxwell's stress tensor, 



TndA, 



(A3) 



where the integration is taken over the surface dA of the 
object and n is the outward normal vector to the closed 
surface S. Maxwell's stress tensor is given by the electric 
and magnetic fields 



E a E p + c 2 B a B p - i(EE + c 2 BB)5 Qf 



(A4) 

Plugging in the expression for the scattered electromag- 
netic field in the above equation, we can determine the 
forces on the dielectric. This method to determine forces 
enables one to calculate trapping of dielectrics and also 
to calculate radiation pressure effects. However, this ap- 
proach cannot be used to determine a full dynamical de- 
scription of the system and its decoherence rates. 



Appendix B: Wigner-Weisskopf with correlations in 
the field 



In this Appendix, an alternative approach to the de- 
scription of the interaction between a single photon and 
a dielectric is given within the Wigner-Weisskopf ansatz. 
In contrast to a direct description with master equa- 
tions, where information about the light fields can be ex- 
tracted via the quantum regression theorem, the Wigner- 
Weisskopf approach directly yields expressions for the 
photonic fields. Analyzing the light emitted from the 
cavity yields information about the mechanical state of 
the system. Following [53] , it is possible to determine the 
occupation of the mechanical mode, and thus, to moni- 
tor the cooling of the system. Complementing the mas- 
ter equation ansatz, in App. |B 1| we solve the equations 
of motion for the coefficients of the density matrix tak- 
ing into account correlations among the free modes of 
the field. In full analogy to the solution of the Heisen- 
berg equations of motion given in Sec. |III A[ the equa- 
tions of motion for the coefficients of the density matrix 
can be demonstrated to be equivalent to the Lippmann- 
Schwinger equation enabling one to use solutions of the 
classical scattering equations. Subsequently, an inhomo- 
geneity in the electromagnetic field is added in App. |B 2| 
and its effect on the scattered fields is investigated. Fi- 
nally, a master equation describing the decay of the cavity 
mode is derived in App. |B2b| 



1. Free photons 

Here, the evolution of a single photon in a plane- wave 
state interacting with a dielectric is discussed. For sim- 
plicity, the motion of the object is neglected for now, 
assuming M — > oo. The coefficients of the wave function 
in the Schrodinger picture are defined by 



|V) = / dkc k (t)al\n), 



(Bl) 



where |fi) denotes the vacuum state. The assumption 
that the object's mass is infinite manifests itself in the 
independence of the wave function of the object's mo- 
mentum state: the effect of the photon's recoil on the 
dielectric is neglected. To obtain the Eqs. of motion, 
we let the homogeneous part of the Hamilton given by 
Eq. ( 10 1, Hb, act on the above wave function. This yields 



ick(t) =WkCk(£) + £ c / dx dk'^k^k' 

Jv{t) J (B2) 

x (c k '(Oe i(k ~ k ' )x + c^(t)e i ( k+k ') x ). 
In order to close this system of equations, we define 

/jM(x,f,i) =i J dVZJk~e- ik(x+f) c k (<), (B3) 
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with h(ic,f,t) = frW(x,r,t) + h^(x,r,t). Sub- 
sequently we multiply both sides of Eq. (B2| by 
i J rfky^tJk exp(— ikx). 



A formal integration over time 
and a transition to the frame rotating at a frequency 
cj , ^(x,?,*)^' = exp(±iw i)^ + H x ^)j yields 



(x, f, t) =h£ ) (x, r,t)+e c [dk [ dx'e^ 

7 Jy(f) 



k(x-x') 



Wk - w + i7 



(B4) 



where ftj^ (x, f, t) is defined in analogy to Eq. ( 16 ). Also 
here, we assume h^ + ' (x, f, t) to be peaked at uj . In or- 
der to solve this differential equation, we have taken the 
slowly-varying approximation, assuming that h(x, f, t)^' 
can be taken out of the integration that is extended to 
t — > oo. A transformation back to the coefficient picture 
thus gives in full analogy to the operators of the electro- 
magnetic field, Sec. Ill A 



Ck(t) 



dk'e 



X <5(k-k') 



7k,k'(^ 



Ck'(O). 



(B5) 



The coefficients contain the information about the full 
dynamical evolution of the system and can be used to 
reconstruct its density matrix. 



2. Cavity field 

In this section the analysis of the previous chapter is 
extended to the more general case, where an inhomogene- 
ity in the electromagnetic field is present. This inhomo- 
geneity is typically a cavity that changes the system's 
mode distribution. Also in this case, the cm degree of 
freedom is treated as a number and its motion in the 
cavity is neglected. We solve the Schrodinger eqs. of 
motion for the coefficients in Sec. |B 2 a| and subsequently 
derive the master equation for the time evolution of the 



cavity mode So in Sec. B 2 b 



a. Solution of the inhomogeneous part 



The wave function including the cavity mode is written 



as 



|V>) = co(t)4|fi) + J dkc k (t)al\Sl} + s \n), (B6) 

where sq is constant and co(t), Ck(i) are time-dependent. 
The Hamiltonian, Eq. Q causing the scattering con- 
sists of two parts, W(f), Eq. (fTTl) describing the coupling 



among the homogeneous modes and -ffesi Eq. (12) de- 



noting the coupling of the inhomogeneity (cavity mode) 



to the free ones. The following equations of motion are 
obtained: 



ick(t) =e c / dx / dk'e^-^V^c^CkKi) 
Jv(f) J 

+ < ,. I dx A /^e- ikX /(x) C0 (i) + CkCk(t) 



v(f) 



2Vn 



and 



V(f) 

■ w c (t), 



ico(t) =e c / dx / dk A /^ e ikx r(x) Ck (t) 



2V n 



(B7) 



(B8) 



where the rotating-wave approximation has been 
taken, which is equivalent to remaining in the one- 
excitation manifold. Using Eq. (B3), we can simplify 



Eqs. (B7), (B8) to obtain after an integration over time 



h$ (x, r, t) = (x, f , t) + d(x, f, t) (B9) 



with 



d(x,f,t) = J dkA(k,f,t)e ikx , 

1 

A(k,f,t) =e c [ dre- i ( w °- Uk)T x 
Jq 



(BIO) 



x / dx', M ^e- ikx '/(x')co(r). 
v(f) V zv o 



The first part of the integration in Eq. (B9 



(BH) 
has been 

carried out under the Markov assumption, which is jus- 
tified as correlations in the electromagnetic field decay 
quickly. No approximation is taken for the time evolu- 
tion of the inhomogeneous part and it is kept in the most 
general form for now. The strategy to find a solution for 
the inhomogeneous case described by Eq. ( B9 1 is to con- 



nect it to the homogeneous one, described in the previous 
section, Sec. IB II 



The solution of the homogeneous case, Eq. ( B4 1 , can 
be formally written in vector-form as 



h = h in + Bh 
1 . 



(B12) 



l-B 



where B describes the scattering operator in matrix form 
and h denotes the continuous vector-representation of 
h(x,r,t). Comparing the homogeneous case to the in- 
homogeneous one, an additional inhomogeneous term is 
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present leading in analogy to Eq. (B12| to 



h = hi, 



h = 



Bh + d 



1-B 



~(h in + d), 



(B13) 



where 1/(1 — B) denotes the solution-operator for the 
plane-wave state. This equivalence facilitates the solu- 
tion of the inhomogeneous system by connecting it to 
the homogeneous one. The system is initially assumed 
to have one photon in the cavity mode and none in the 
homogeneous modes, c k (0) = 0, such that hj n = and 
Eq. (B9) can be solved as 



These results are in accordance with Fermi's Platinum 
Rule, the extension of Fermi's Golden Rule to all orders 
in multiple-scattering processes. 



From the coefficients to the master equation 



Starting from the wave function Eq. (B6 1, the system's 



density matrix is obtained by tracing out the environ- 
ment and can be written as 160 



Ps(t) 



/i (+) (x,f,i) =d ta (x,f,t) + . 

w k 



dk I dx'e- ik(x - x ' } 

V(f) 



V s oCo(*)* 
Taking the derivative gives 



w k - w + 17 



d(x,f,t), 



(B14) 



Ps(t) 



s c (ty 



s^c (t) 
1-Mi)| 2 



di 



\coW 



where d in (x,f,t) is defined in analogy to Eq. dBb. Tak- Usin S Ec t s - d B18 h d B19 D finall y y ields 



ing the inverse transformation yields the solution for the 
coefficients c k (t) 



Ps 



A L a ao,Ps 



c k ft) = J dk'e-'^x 
x U(k-k'H 



with 



T k ,k'(f) 
Wf - w + i7 



£ L [ps] = « feaops&l - [ala ,ps)]j 



(B20) 



(B21) 



(B22) 



(B23) 



(B15) 



Plugging Eq. (B14| back into Eq. (B8), we can close the 



equations of motion. After taking the Markov approxi- 
mation assuming that the system does not change signifi- 
cantly during the interaction with the environment, such 
that Co(f — t) w co(t). Using some standard relations for 
the scattering operators, the time evolution of the inho- 
mogeneous mode in terms of transition operators is given 
by 



where A L denotes a shift of the energy levels, and k de- 
scribes the decay rate of the cavity photons due to losses 
effected by the presence of the dielectric. This master 
equation and its decay rates are equivalent to the result 
for the light fields obtained in Sec. |III| 



co ft) 
c ft) 



= - / dk\T Kc (r)\ z 



7r5(w k — wo) + iT 5 - 



1 



w k 



- w 
(B16) 

The effect on the light field can be determined approxi- 
mating Tk c (f) ss 7^,0(0) thus neglecting the effect on the 
cm mode. Integration gives 



c (i) = exp((-/c + iA L )f)c o (0) 
with the decay rate 



and a Lamb shift 



= -5R 



L c> 



A = 



co ft) 
c ft) 



CO ft) 

co ft) 



(B17) 



(B18) 



(B19) 



Appendix C: Contributions to the master equation 
in the Lamb-Dicke regime 

In this Appendix, we demonstrate in some detail which 
terms of the master equations discussed in Sec. |III C| can 
be neglected within the Lamb-Dicke regime under strong 
driving. We make in particular use of the analytic prop- 
erty of the transition matrix elements and their deriva- 
tives. As an example, we discuss the contribution cx f 
given by 



H 



shift 



a 2 / dk(^(w k -wo)S[r^ c (r)|, =0 7^ k (0)] 



K[72, c (f)|, =0 7^ k (0)])f. 



(CI) 



V- 



w k - w 



Under the assumption that t(w k ) = 7^ c (f)\f=oT* k (0) is 
analytic in w k , the Hilbert transformation can be used to 
show that 



1 

w k - w 



-»[t(wk)] = -9f[t(w )] 



(C2) 
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and consequently H ] 



shift 



0. This transformation is also 



used to simplify the renormalization contribution oc r in 

sec. mrei 

Another argument holds when neglecting the incoher- 
ent contribution arising from terms oc aoa*f. This part 
describes decoherence of the mechanical and the light de- 
gree of freedom 



£S[ ps ] =2T g (a p s (b + V) - [a (b + tf),p s ]) 
.;((b + fr)p s al-[al(b + P),p s } 



■ g\ 
■2T*„ 



(C3) 



x a* / dk6(u, k - ^o)T c : k (0)r k '(f)|, =0 . (C4) 



where 



This contribution can in general be neglected for the cm 
degree of freedom, as it is suppressed by 1/a compared 
to Eq. (38 1. Requiring that at] <C 1, we can also neglect 
the effect of this contribution on the cavity mode. This 
requirement becomes clear comparing Eq. ( C3 ) to the 



decay of the cavity described by Eq. (41 ) 



Appendix D: The small particle limit 

In this section, we consider the important limit of the 
general theory, where the dielectric is smaller than the 
cavity wavelength. In this case, it is justified to ne- 
glect couplings among different bath modes as described 
by Eq. ( |TT| ). While this derivation has been analyzed 
in [331 161) . we demonstrate in this Appendix that the 
same master equation can be obtained starting from the 
general description, Eq. (26 1, and considering only terms 
of the lowest order in the Born series of scattering theory. 
Let us now discuss the different contributions and how 
they are modified: 

1. All terms in H$ remain the same, the trapping fre- 
quency Cb t can be determined from Eq. (J9j) , which 
for objects much smaller than the wavelength is 
given by 



(Dl) 



where p is the material's density, / the laser inten- 



sity of the optical tweezer and Wt its waist. Also 
the optomechanical coupling g is given by Eq. ^ 
simplifying to 



-X - 



4cVv ' 



(D2) 



where V is the object's volume. Note that the ex- 
pressions for Hs are not affected by neglecting the 
coupling among the bath modes as these quantities 
are determined only by system operators. 
The renormalization of the optomechanical cou- 
pling and the trapping frequency are obtained by 
considering only the lowest order of the Born series 
7k c (f) ~ Vkc(f) i n the expressions for A and g ln 
given by Eqs. p7|), Q. This gives 
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° = -e c k 2 R 2 c 



(D3) 



with R being the sphere's radius. The renormal- 
ization of the trapping frequency is obtained by in- 
serting the trapping mode and leads to 



A M <° - -e c k 2 R 2 ul 



(D4) 



3. The same procedure, namely taking the lowest or- 
der in the Born series by setting T c ,k(r) w V c ,k(^) 
is also applied to obtain the decoherence rates. For 
the cavity decay rate this gives 



K o _ e c K v c 
24ttV c 

and for the recoil heating of the cm we obtain 
Pt . n ph c 



(D5) 



pO _ € c K V 



Qirpuj t \uj L nW 2 2V C 



(D6) 



where the first term describes decoherence due to recoil 
heating by photons from the tweezers and the second 
term by photons of the cavity mode. Comparison to |23j 
shows that this result is in agreement with the what is ob- 
tained when directly taking the Born Markov approach, 
neglecting interactions between bath modes and deriv- 
ing the master equation in the standard way. (Note that 
there are some differences in prefactors as we take into 
account the polarization of the light here.) 
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